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Overview of Research 
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Targeted Platform 


JPL Radio 

¢ 66-MHz SPARC processor 
¢ 2.x Xilinx Virtex-Il FPGAs 
¢ S-band transceiver 

¢ L-Band (GPS) receiver 


Picture Source: Reinhart, R., Wireless Innovation Forum 
Technical Conference, 2013 
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Interacting Multiple Model (IMM) Background 
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Picture Source: Sathyan, T. and Kirubarajan, T., /EEE Transactions on Circuits and Systems, 2006 
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Interacting Multiple Model (IMM) Background 


Step Step Iterations | Expression Executed 

Number Type Needed 

MIX01 SF 1 fi(k + 1|k) = pw" (k) x at 

MIX02 SF 1 a xdiag((k)) 

MIX03)— SF 1 u(k) = (w" xdiag(u)(k))” x (diag(f(k + 1|k)))~* 

MIX04 SF 1 Xo(k|k) = X (klk) x p(k) 

MIX05 DF R? X jo(k|k) — Xi (k|k) 

MIX06 DF R? (Xjo(k|k) — Xi(k|k)) x (Xjo(k|k) — X4(k|k))? + Pi (klk) 

MIX07 DF FR Pyo(klk) = [(Xjo(Klk) — Xa(k|k)) x (Xjo(klk) — XG (k|k))? + Pl 
x diag (1:3 )(k) + G 

KF0l DF R X;(k+1|k) =F x Xjo(k|k) 

KF0O2 DF R F x Pyo(k\k) 

KF03 DF R P;(k + 1|k) = F x Pyo(k|k) x FT +Q 

KF04 DF R P;(k +1|k) x HT 

KF05 DF R S;(k+1) =H x P(k+1|k) x HT +R; 

KF06 DF R K(k +1) =P,(k+1|k) x H™ x (S;(k+1))7! 

KF0O7 DF R P;(k + 1k +1) = Pi(k + 1|k) — Kj(k +1) x (Pj(k + 1k) x HP)? 

KF0os DF R e;(k +1) =2(k+1)—H x X;(k+ 1k) 

KF0O9 DF R Xj(k + 1)k +1) = Xj(k+ 1k) + Kj(k 41) x e(k +1) 

COM01 DF R el (k+1)x S>'(k+1) xe (k+1) 

COM02 SC R —4eP(k+1) x S>'(k+1) x ej(k +1) 

COM03 SC R exp(—se; (k +1) x Si" (k +1) x e;(K+1)) 

COM04 SC R |S;(k+1)| 

COM05 SC R (|S;(k + 1)|)72 

COM06 SC R L ,(k +1) = (2m)? (|S;(k + 1)|)~2exp(—ged (k +1) x S51 (k +1) 
x ej(k+1)) 

COM07 SF 1 diag(L(k +1)) x ji’ (k + 1k) 

COM08 SF 1 ji(k + 1|k) x L7(k+1) 

COM09 SF 1 u(k +1) = (diag(ji(k + 1|k) x L’ (k +1)))~! x (diag(L(k + 1)) 
x fi" (k + 1|k)) 

COM10 SF 1 X(k+1)k +1) =X(k+ 1k 41) x p(k 41) 
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Implementation Trade Studies 
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Measurement (Noisy) 


Prediction (Estimate) 


Picture Source: Faragher, R., /EEE Signal Processing Magazine, 2012 
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Schur Complement Architecture 
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processing element 
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out element 


processing element 


Slave 1 processing element controller 


processing element 


Slave K processing element controller 


Sudarsanam, A. et al., JET Computers & Digital Techniques, 2010 
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Schur Complement Analysis 


Multiplier Addition Divider Cycles Divider | Occupied Slice Flip 4-Input Max Clock 
Latency Latency Per Operation Latency Slices Flops LUTs_ | Frequency (MHz) 
1 1 14 14 | 541 1299 00.742 
3 i) 14 14 S26 925 1340 (BS ss: 
a 7 26 ra 1090 1328 152.879 
6 7 26 27 1201 1343 175.389 
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IMM Filter Architecture 
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IMM Analysis 
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imm_filter 
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Example Scenario 


State Model 
Constant Acceleration Model 


Available Measurements 
SNR, SIR 


System Models 
Clear Sky, Rain Fade, 


lonospheric Scintillation, 
Tropospheric Scintillation, 
Multipath Fading 
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Filter Rate (kHz) 
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IMM Analysis 
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10 15 
Number of Processing Elements (PEs) 


—+— KFs = 2, X Elements = 2, Z Elements = 3 
—+— KFs = 2, X Elements = 2, Z Elements = 4 
—— KFs = 2, X Elements = 4, Z Elements = 2 
—+— KFs = 6, X Elements = 2, Z Elements = 3 
A,B,C = PE Controllers, 
PEs Per KF Controller, 
PEs Per IMM Controller 
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Summary 
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If youre interested, contact me at: 


Tim Hackett: tmh5344@psu.edu 
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